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Про один підхід до ідентифікації особи 
за контуром профілю носа 


У статті розглядається підхід для ідентифікації особи за контуром профілю носа. Для знаходження 
простору ознак пропонується використовувати апроксимацію за допомогою В-сплайнів. Алгоритмічна 
реалізація підходу використовує моделювання структури бази даних у вигляді п- вимірного куба. 


Вступ 


Біометрична ідентифікація побудована на порівнянні деяких параметрів (як фізіо- 
логічних, так і психологічних), що властиві людині та занесені до бази даних. У випадку 
фізіології для аналізу та порівняння використовуються такі постійні та незмінні параметри 
людини, як фізичні характеристики: риси і форма обличчя, особливості очей, рук - долоні 
та пальці (папілярний узор, топографія, форма), венозний рисунок, теплова картина тощо. 

Останнім часом у дослідників викликає інтерес і вважається перспективною іденти- 
фікація людини за формою носа (1 |. Одним із варіантів застосування такого способу вста- 
новлення особи є виявлення злочинців, знятих на камери прихованого спостереження. 

Алгоритми локалізації облич на зображенні та виділення контурів на обличчі до- 
сить розвинені. Метою даної роботи є задача ідентифікації особи за контуром профілю 
носа, за умови, що такий контур вже отримано. Матимемо на увазі, що процедури отри- 
мання контуру (навіть за умови його віднормованості) не гарантують повний збіг в од- 
нієї і тієї ж особи. Це випливає із особливостей як зображення, так 1 цих процедур. Тому 
розв'язок задачі ідентифікації має це враховувати 1 вирішувати. 


Формалізація підходу 
Алгоритми обробки зображень дозволяють отримати контур носа як множину 
точок: 
2 «(2 -(оуфієо-М -1ху, є ВІ. а) 
Подальша робота з такою множиною точок можлива за умови проведення 
відповідного нормуючого перетворення: 


Уа В (2) 
де Фа и, зберуріе 0--М -1х, у, є ом. 
Пропонується (2) отримати за допомогою: 
д,-з24, (3) 
де А - матриця (3 х3) перетворення, яка визначається як добуток матриць переносу, 
повороту та масштабування; а, - (х, уві) 2 (х, уд 1) - точки в однорідній системі ко- 


ординат. Отримаємо матрицю А в аналітичному вигляді. Задача полягає у визначенні 
матриць: переносу (Т) - для переносу кривої-контуру в початок системи координат, 
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повороту (КВ) - для суміщення початкової та кінцевої точок кривої з оссю Х та 
масштабування (5) - для розподілу точок кривої між нулем та одиницею: 


І 0 0 сок(а) -зіп(а) 0 5. 0 0 
т-| 0 1 0|,В- | віл(а) сока) 0|,5-10 з 0, (4) 
не ліяфь м 0 0 1 0 0 1 


де сова - (ха ві о; зіп а «(уз є 0; О- ЯІРЯ зда зу, -Уйа у - від- 


стань між крайніми точками кривої /; 5 - у 0,350, зха, 7 Хо: 
Матриця перетворення А матиме вигляд: 
5 со5 (0) -58іп (0) 0 
АЗеТ-В-5- 58іп (а) 5 со5(аг) че 
-5( в віє) -- ху соз(аг)) | -5 (о 8іп(а)- уд со5(а)) 1 
, 5 

уо -0, 0,0) 0 
а 0,0, 0) уо 0 


- (о 0, ку 02/00, 0) (хо :0, - у, 0,)/(0,:0) 1 
де 0,зУд 17 Ус: 

Процедури визначення точкових кривих, що відтворюють контури носів різних 
людей, продукують різну кількість точок ( М ). Подальші перетворення можливі за 
умови рівності кількості точок ( М ) у кривих-контурах. Досягнути рівності кількості 
точок можливо, скориставшись кусково-лінійною інтерполяцією: 

хахаівохоИМОЇ 


ічі1 


- т, |З ту Ж лоз ра гео ЗЕ 2 
я ом зх а 9 ЗДО 


(6) 


де 2232, 17)» із1,5М-1, є), : -); х, У - вхідна (неінтерпольована) множина. 
Поставлена мета полягає у знаходженні відображення множини точок (Г)) у 
множину ознак (Р): 

р-Р, (7) 
де Ро і рт збо У; і 2р8М У , М «« М , елементи якої єдиним чином визначають контур 
носа. Із фізіогноміки, живопису |2| та криміналістичної експертизи |3| відомо, що про- 
філь носа має декілька зон-особливостей, комбінація яких і дає відомі контури носа 
(грецький, римський, орлиний тощо). До таких зон відносяться основа носа, кінчик 
носа, спинка носа та ще можна віднести перенісся. Перенісся враховується для більш 
повної типізації носів. Топологія кривої, яка відтворює контур носа, має яскраво вира- 
жені чотири зони особливостей. Звідки, для перетворення (7)- М - 4. 

Перетворення (7) пропонується знайти, використовуючи В-сплайн-апрокси- 
мацію. Задача В-сплайн-апроксимації є задачею підгонки В-сплайн-кривої з К-опор- 
Т -- -. уд 
ними точками ра Нр до точкової кривої а з (4, ; сао) , де М»К (зазви- 
чай М2»К) для значень параметра мо... |. Така задача апроксимації приводить 
до перевизначеної системи лінійних рівнянь: 


Ми) раза, (8) 
де М, (п) - В-сплайн-базисна функція. 
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Підхід, вибраний для формування В-сплайнів, полягає у тому, щоб визначати 
сплайн в термінах множин базисних функцій, кожна з яких відмінна від нуля тільки 
на інтервалі в декілька вузлів (значень параметра и). Кожна функція М; (ш) є полі- 
номом заданого степеня на інтервалі в декілька вузлів та дорівнює нулю за межами 
інтервалу. Множина функцій їм (ш)) утворює базис для заданої послідовності пара- 


Ії 


метрів и і заданої степені полінома. 
Існує багато способів визначення базисних функцій ( М, (ш)), але особливе місце 


належить рекурсивним функціям Кокса-де Бура (41: 


І, якщом, хи хи, 


Хот . 
і 0- у іншому випадку 
і (9) 
и-и иго-и 
а К ка 
Ма З Муачи) В Міна (и) 
Чеа" Чу Чан З Мр 
де й, Цій, - пОСЛІідОВНність вузлів, така що: 
й ЗНОСИН оЗИ, зи охо 


Кожна функція із першої множини, М, , , постійна на одному інтервалі та дорів- 
нює нулю за його межами; кожна функція із другої множини, М, ,, лінійна на двох 
інтервалах та дорівнює нулю за їх межами; кожна функція із третьої множини, М, », 


має вигляд квадратичної кривої на трьох інтервалах та дорівнює нулю за їх межами. 
В загальному випадку функція їх множини М, , має відмінне від нуля значення на 


а інтервалах між и, і и, ,,, 1 має вигляд полінома степені а на кожному із цих 


інтервалів, причому у вузлах забезпечується параметрична неперервність класа С". 
Сплайн-крива, утворена як зважена сума таких базисних функцій, буде лежати все- 
редині опуклої багатокутної оболонки, так як для базисних функцій виконується 


мМ-1 
умова: У Ми) -1,0х« Ми) «1. Але від того, що кожна функція М, , відмінна від 
із0 


нуля тільки на а--1І інтервалах, то кожна опорна точка впливає тільки на ту діль- 
ницю сумарної кривої, яка лежить всередині оболонки, утвореної 4-1 опорними 
точками. 

Множина базисних функцій В-сплайна визначається заданим степенем сплайна та 
вектором вузлів. Не накладено ніяких обмежень на значення у вузлах, крім як и, Хи, |. 


Якщо вважати, що результат ділення на ноль при обчисленні по рекурсивним формулам (9) 
дорівнює 1, то можна мати кратні вузли (ті, що повторюються). При повторенні вузлів 
(кратності вузлів) проявляється ефект наближення сформованого В-сплайна до опорної 
точки, яка асоціюється з цими вузлами. Якщо кратність вузла буде дорівнювати а --1, 
то В-сплайн степені а буде проходити через відповідну опорну точку. 

Формування кубічного сплайна проводиться з використанням наступного масиву 
вузлів: 10, 0, 0, 0...,1,1,1, й . Наслідком цього є належність першої та останньої опор- 


ної точки самій кривій. Тому для отримання потрібних для аналізу п точок потрібно взяти 
п-2 точки. 

Базисні сплайни заданої степені є лінійно незалежними функціями та утворюють 
базис у функціональних просторах, що можна використовувати для відображення з їх 
допомогою інших функцій цих же просторів. Виходячи з цього, для отримання набли- 
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женого розв'язку перевизначеної системи лінійних рівнянь (8) можна скористатися мето- 
дом найменших квадратів (51, 161: 
Т Т 
М М'рем' 4, (10) 
звідки 
-1 
ре |м'м)"-м'а. (11) 
Задача відтворення за визначеними за допомогою (11) опорними точками вхідної 
кривої для профілю носа потребує від 9 до 16 точок. Це пов'язано з тим, що, крім озна- 
чених вище чотирьох зон-особливостей, у профілів носів є ще ряд незначних (з точки 
зору особливостей) кривизн. Пропонується визначати чотири опорні точки. З точки зору 
апроксимації, точного відтворення вхідного профілю не отримаємо, але у цих чотирьох 
точках будуть сконцентровані особливості саме в потрібних нам чотирьох зонах. 
Застосування перетворення (7) для множини контурів носів, які треба ідентифікувати, 
приводить до отримання чотирьох множин параметрів для зон особливостей профі- 
лю носа: 


Р! - Ку/, та), (2) 
де і) - 1,... 4 - для чотирьох зон; 1- ,..., І,, де І, - кількість профілів носів різних 
людей; Іа - номер, який однозначно ідентифікує людину. 

При формуванні множин (12) використовуються тільки у-координати точок-ознак. 
Це випливає з того, що перед перетворенням (7) множина точок (контур носа) піддалася 
кусково-лінійній інтерполяції (6). Як наслідок цього, х-координати точок-ознак різних 
контурів будуть однакові, тому їх можна виключити з подальшого розгляду. 


Тоді задача ідентифікації полягатиме у пошуку такого ла є реа нон, 
де Бі КУ та ув - Р|«А,і-1-2 є Р, 


де Б- є Ка 41 - множина ознак вхідного контуру для пошуку. 


Радіус пошуку В підбирається із врахуванням деяких незбігів між вхідним конту- 
ром та контуром, що записано раніше. 


Алгоритмічна реалізація підходу 


Основна ідея підходу полягає у відображенні множини (12) у координатний 
простір |7|. Як обчислювальну платформу пропонується використати 501,-сервер. 
Використання традиційної реляційної моделі з властивою їй декомпозицією досить 
громіздке. Беручи до уваги, що множини (12) досить сталі, для проектування бази 
даних пропонується структура предметної області у вигляді п-вимірного куба |71. 


А 13) мантіх: таблица 
А|В| Оаїа! | Оаїа? | 

ДЦ | 0 1 Прізвище 
| | 0 2їмя 
| | 0 3 Побатькові 
| | 0 4 Датанародження 
| | 0 11 Точка-ознака 1 (у-координата) 
| | 0 21 Точка-ознака 2 (у-координата) 
| | 0 31 Точка-ознака 3 (у-координата) 

і | | 0 41 Точка-ознака 4 (у-координата) 
| | 0 12 Максимальна відстань між точками у множині 1 
| | 0 22 Максимальна відстань між точками у множині 2 

В | | 0 32 Максимальна відстань між точками у множині 3 
| | 0 42 Максимальна відстань між точками у множині 4 
Рисунок І - 2-вимірна структура Рисунок 2 - Довідники 
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Предметна область задачі має єдиний об'єкт «людину» з рядом параметрів. Пріз- 
вище, ім'я, дата народження, координати точок-ознак -- параметри об'єкта «людина». 
Всі параметри унікальні, тому не потрібно реалізовувати підлеглих структур. Для од- 
ного об'єкта з декількома параметрами п-вимірна структура вироджується у 2-вимірну 3 
вимірами А та В (рис. 1). Розглянемо більш детально виміри. При А - 0 значення даних 
на осі В є, по суті, довідниками (рис. 2). Значення виміру А від / до І, (де І, - кількість 
людей) призначені для фіксації інформації по конкретних людях (рис. 3). Для збере- 
ження, модифікації та пошуку даних (4 ) використовуватимемо: 


4; - (а,,Б)), (13) 
де а,є То І , Б,є П, 2,3,4,11,21,31, 41 , 4, - значення, які знаходяться за ко- 


ординатами (а,,5,). 


3 Майтіх: таблица 
А Оаїа! Оаїа? 
3 маітіх: таблица | РА|В | Оаат | Оаа? | 
|У 1 1 Їванов 
перу оре! Оаїаї | | 1 2іван 
1 Ї|ванов | | 1 З Іванович 
1 Т 1. 4101.03.1963 
2 1 Петров Їм ери 
3 1 Сидоров | ШОЕ 0491 
ЖЕ 0.05 
| 1 4 40.03 
Рисунок 3 - Вимір А - інформація Рисунок 4 - Інформація 
по людях по одній людині 


Можливі комбінації а, і Б, означатимуть: (а, 1) -- прізвище, (па, 2) - їм'я, (па, 3) - 
по батькові, (та 5 4) - дата народження, (та ,1 1) - у значення першої точки-ознаки, 
(па ; 21) - у значення другої точки-ознаки, (та ,3 1) - у "значення третьої точки-о0з- 
наки, (та , 41) - у-значення четвертої точки-ознаки ( па є ет 1 , де І, - кількість 


людей. Беручи до уваги, що запропонована структура реалізовуватиметься за допо- 
могою реляційної СУБД, і пошук виконуватиметься засобами цієї ж СУБД, то для опти- 
мізації варто область значень подати двома стовпчиками: Дата! -- для нечислової інфор- 
мації (текст, дата) та Да!а2 - для числової інформації. Тоді стовпчик Дагад можливо 
проіндексувати, що прискорить пошук, організований засобами реляційної СУБД. 
Приклад заповнення інформації по одній людині наведений на рис. 4. 

Відображення потрібної нам предметної області потребує лише (на відміну 
від реляційної моделі) однієї таблиці. Засоби роботи з цією таблицею (58іогед- 
процедури з 501,запитами) також суттєво спрощуються. Запис та модифікацію 
довільного рядка таблиці реалізує наступна 5іогед-процедура (лістинг 1): 

Лістинг 1 - Зіогед-процедура для запису та модифікації рядка у таблиці 


СЕКАТЕ РКОСЕРОВЕ аро.Іпзегі 

ад іпе,Фв іпе,Фрабаї Уахгспаг (50),ФРага2 Еісаї 
АЗ5 

Ів 


(5ЕШЕСТ СООМТ(х) ЕБОМ Маєнгіх МНЕВБЕ (А-ФВА)АМО (В-08В))»0 
ВЕСІМ 
ОРРАТЕ Маєгіх З5ЕТ А-ВА,В-ЄВ,Драса1-вРасга1,раса2-в8Раса2 ПНЕВБЕ (А-8ФА) АМО (В-8В) 
сиро 
ЕТЗЕ 
ВЕСІМ 
ІМ5ЕВТ ІМТО Маєгіх (А,В,рака1,рага2) МАЦОК5 (Ф8А, «В, Єрака1, "Раба?2) 
сиро 
со 
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Визначення особи за чотирма вхідними точками-ознаками профілю потрібує до- 
датково розрахунку максимальної відстані (радіусу пошуку) між двома сусідніми точ- 
ками для кожної з чотирьох множин. Розрахунок реалізовано за наступним алгоритмом: 

- вхідна інформація: значення у-координати нової точки і-ї множини у, 


(і р 1,:,4); 
- шукаємо першу точку (у,) за умовою у, » у, та розраховуємо відстань: 
а «у, 
- шукаємо другу точку (у,) за умовою у, «у, та розраховуємо відстань 
4, З» У, 


- вибираємо максимальне значення відстані 


2 


2 


-- тах(а,, 4, ); 

- якщо знайдене значення максимальної відстані більше за записане раніше 
(для відповідної множини), то модифікуємо його. 

Наведений алгоритм реалізує відповідна 5іогед-процедура (лістинг 2). Вона 
виконується після кожного запису або модифікації значень точок-ознак. Її результа- 
том є модифікація рядків з координатами (0), 12), (0,22), (0,32), (0,42), поле Фага? яких 
містить поточне значення максимальної відстані для кожної множини. 


тах 


Лістинг 2 - Зіогед-процедура для визначення максимальної відстані 
між значеннями у відповідній множині ((2)В) точок-ознак 


СБЕАТЕ РЕОСЕРОВБЕ аро.МодіїісаніопРрівнкапсе 
ав іпо,ФУп Єбїоає, ФРака2 їїсаї 

А5 
ІЕ (ЗЕГЕСТ СООМТ(х) ЕКОМ Маїргіх МНЕВЕ (А»0)АМРр (В-ЙФВ))21 
ВЕСІМ 

ОЕСІЬАВЕ Шар їбіоаєг, ФМах їісає 
ЗЕФЕСТ Фре (5ЕТЕСТ МАХ (у) ЕБОМ ( 
ЗЕФТЕСТ ТОР 1 АВ5 (Драга2-вб8Уп) у ККОМ Маєгіх Мпеге А-0 АМР В-8В АМО Раба2»вуп 
ОМІОМ 
ЗЕФТЕСТ ТОР 1 АВ5 (рага2-вб8Уп) у ККОМ Маєгіх Мпеге А-0 АМР В-8В АМО Раба2«ауп)) 
ЗЕТГЕСТ ОМАХ-(ЗЕТЕСТ РАТА2 ЕБОМ Мабгіх ЙНЕБЕ В-й8В АМР А-б0) 

ІК аМах«Ййр 
ВЕСІМ 
ОЮРРАТЕ Макгіх ЗЕТ Расга2-0йРр ЕБОМ Мабгіх ТНЕБЕ (В-ЙФВ) АМО (А-0)) 
ЕМО 

ЕМО 
со 


ші іє 


Ідентифікація особи за точками-ознаками профілю носа виконується за наступ- 
ним алгоритмом: 

- вхідна інформація: чотири у-координати точок-ознак шуканого контуру; 

- для кожного значення вхідної у-координати у відповідній множині шукають- 
ся значення, що попадають в окіл (1/- максимальна відстань між двома значеннями у 
множині); 

- шукається перетин отриманих множин за умови наявності у них значень однієї 
особи; 

- якщо такий перетин не пустий - шукане знайдене, в іншому випадку - шукане 
не знайдене. 

Наведений у лістингу 3 код 5іогед-процедури реалізує ідентифікацію особи за 
описаним алгоритмом. 
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Про один підхід до ідентифікації особи за контуром профілю носа 


Лістинг 3 - Зіогед-процедура для ідентифікації особи 


4Б 


СЕЕАТЕ РЕОСЕРОВК аро.Ідепніїісабіоп 

-- вхідні значення 

ау1 Еїоаїк, ФУ2 їіоає, ФУ3 Еїоаг, ФУ4 Еїсаєм, 

-- вихідне значення 

саІпа іпіь 
АЗ5 

-- визначаємо максимальні відстанні між значеннями для 4 множин у-координат точок-ознак 
ОЕСІЬАБЕ ФРріз5і1 їіоаїг, ФРіз5с2 їїоаїк, Ррі5і3 їїоаї, ФРрі5і4 їїсаї 

ЗЕФТЕСТ ФРрі5:1- (5ЕТЕСТ Ракба2 ЕКОМ Маїкгіх МНЕБЕ А-0 АМОЬ В-12) 

ЗЕФЕСТ ФРрі5:2-(З5ЕФЕСТ Равка2 ККОМ Маїкгіх "МНЕБЕ А-0 АМЬ В-22) 

ЗЕФЕСТ ФРріз5ь3-(З5ЕТЕСТ Рака2 ЕКОМ Маїкгіх МНЕБЕ А-0 АМР В-32) 

ЗЕФЕСТ ФРрі5ір4-(5ЕТЕСТ Рага2 ЕКОМ Маїгіх МНЕБЕ А-0 АМЬ В-42) 
--визначаємо 4 підмножини з кандидатами на ідентифікацію 

дагор каріе йТ 

СЕКАТЕ ТАВІЕ йТІ (|А іпео| МОТ МО11,) 

МЗ5ЕВТ ІМТО ЙТІ 5ЕГЕСТ А ККОМ Маїбгіх Мреге АВ5 (Ррака2-У1)«-Різі1 АМР В-1 
дагор Ркаріе ЙТ2 

СВЕАТЕ ТАВІКЕК йТ2 (|А іпе| МОТ МОГ) 

МЗ5ЕВТ ІМТО йТ2 5ЕГЕСТ А ККОМ Маїгіх Мреге АВ5 (Рака2-Уу2)«-Ррізі2 АМР В-2 
дгор Ркаріе йТ3 

СВЕАТЕ ТАВІКЕ йТ3 (|А іпе| МОТ МО11) 

М5ЕБТ ІМТО ЙТ3 5ЕГЕСТ А ККОМ Маїкгіх Мреге АВ5 (Рафка2-Уу3)«-Різь3 АМР Ве-з 
дагор каріе їТ4 

СЕКАТЕ ТАВІЕ йТ4 (|А іпео| КОТ МОГ) 

М5ЕВБТ ІМТО йТ4 5ЕГЕСТ А ККОМ Маїкгіх Мреге АВ5(Ррага2-.4)«-Різн4 АМрР В-д 
-- визначаємо результат ідентифікації як перетин 

ЗЕТЕСТ аФІпда- (ЗЕТЕСТ йЙТ1.А аз Іпа КВОМ 8Т1,172,1Т3,йТ4 ЙНЕБЕ Іпа-йТ2.А АМО Іпда-йТ3.А АМР 
па-йТА4.дА) 

со 

Експеримент 


Для перевірки спроможності запропонованого підходу була побудована множина 


із 17 контурів профілів носів (рис. 5). 


Рисунок 5 - Контури профілів носів 


Для побудови цієї множини фотографічне зображення профілю людини переводи- 
лось до контурного вигляду, скелетизувалось, отримувались точки контуру носа та за допо- 
могою перетворень (5) 1 (6) отримувались контури носів для подальшого аналізу (рис. б). 


Рисунок 6 - Отримання контуру носа 
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Застосувавши перетворення (7) для цих контурів носів, було отримано чотири 
множини (12). Рис. 7 містить графічне відображення цих множин. 


05 


0.4 


оз 
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а Спинка 


ба - Перенісся 
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о 


-о.1 
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Рисунок 7 - Множини параметрів-особливостей для 17 носів 


Висновки 


Запропонований підхід показав свою ефективність для вирішення поставленої 
задачі. Подальші дослідження будуть направлені на вдосконалення підходу: тесту- 
вання на більшій множині, яка б була репрезентативною для даної території; розгляд 
контурів носів близьких родичів тощо. 
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А.В. Бармак, Ю.В. Крак, Ю.В. Кривонос 

Об одном подходе к идентификации личности по контуру профиля носа 

В статье рассматриваєтся подход для идентификации личности по контуру носа. Для нахождения про- 
странства признаков предлагаєтся использовать аппроксимацию с помощью В-сплайнов. Алгоритмическая 
реализация подхода используєт моделирование структурьі базьт данньгх в виде п-мерного куба. 


О. У. Вагтак, и. У. Ктак, Ти. С. Ктууопо5 

Абоиі опе Арргоасі ої Нитап Ідепіййсабіоп Бу пе Хо5е Сопідиг Ргоїіе 

Ап арргоасії ої питап ідепіїйсаноп Бу Ще позе сопіоиг ргобіє 15 гезеагсред іп Ше рарег. В-5ріїе5 арргохітайоп 
15 виреезісд Гог декегтіпіпє Пе Кеаїите 5расе. АЇдогійфтіс шпріетепіайоп ої їбе арргоасі ц5е5 дакабазе 
зігистиге їп ре Бога ої п-Фітепзіопа! сибе. 


Стаття надійшла до редакції 19.07.2010. 
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